Nonlinear nonequilibrium quasiparticle relaxation in Josephson junctions 
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Nonequilibrium electrons in superconductors relax and eventually recombine into Cooper pairs. 
Relaxation is facilitated by electron-boson interaction and is accompanied by emission of nonequi- 
librium bosons. Here I solve numerically a full set of nonlinear kinetic balance equations for stacked 
Josephson junctions, which allows analysis of strongly nonequilibrium phenomena. It is shown that 
nonlinearity becomes significant already at very small disequilibrium. The following new, essentially 
nonlinear effects are obtained: (i) At even-gap voltages V = 2nA/e (n — 2,3,...) nonequilibrium 
bosonic bands overlap. This leads to enhanced emission of Q = 2A bosons and to appearance of 
dips in tunnel conductance, (ii) A new type of radiative solution is found at strong disequilibrium. 
It is characterized by the fast stimulated relaxation of nonequilibrium quasiparticles. A stack in 
this state behaves as a light emitting diode and directly converts electric power to boson emission, 
without utilization of the ac- Josephson effect. This leads to very high radiation efficacy and to 
significant radiative cooling of the stack. The phenomenon can be used for realization of a new type 
of superconducting cascade laser in the THz frequency range. 

PACS numbers: 74.40.+k 74.50.+r 42.55.Px 85.60.Jb 78.45. +h 



Nonequilibrium phenomena in Josephson junctions 
(JJs) are central for operation of many superconduct- 
ing devices 0, [3, HJ and are important for understand- 
ing fundamental properties of superconductors [J, H, [3] . 
Nonequilibrium quasiparticles (QPs) decay by emitting 
bosons, with which they are interacting. For example, 
predominantly phonons are emitted in conventional low- 
T c superconductors (LTSC) with electron-phonon cou- 
pling mechanism 0, H, H, 113]. Nonequilibrium boson 
emission was also observed in high-T c superconductors 
(HTSC) [HG3, although coupling mechanism in HTSC 
is not yet confidently known. Analysis of the content of 
such emission provides a possibility for direct determina- 
tion of the coupling mechanism in HTSC [13] ■ 

Nonequilibrium effects in stacked intrinsic Josephson 
junctions (I JJs), naturally formed in layered HTSC such 
as P^S^CaC^Os+a (Bi-2212), can be much larger than 
in LTSC JJs |l3[ . This is caused (i) by the atomic thick- 
ness of Cu02 planes with a very small total density of 
states (DoS); and by stacking of junctions which (ii) pre- 
vents QP leakage from electrodes and (iii) leads to cas- 
cade amplification of boson population upon sequential 
QP tunneling [13] • The two latter factors are well known 
in semiconducting heterostructures and lead to a 10 4 -fold 
decrease of the threshold current density in super/lattice 
lasers, compared to single p-n junction lasers [14|]. I ar- 
gue that a similar dramatic enhancement of nonequilib- 
rium effects could occur in stacked IJJs and may facilitate 
population inversion required for lasing [l3t [ . The IJ J 
laser could provide a tunable emission in the frequency 
range ~ 0.3-7 THz and is considered as one of the candi- 
dates for filling the "THz gap" between microwave and 
infrared technologies [13, [13] • 

Understanding LTSC device performance at low T and 
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boson emission from IJJs requires analysis of strongly 
nonequilibrium states. The appropriate microscopic for- 
malism was developed in the BCS work [17]]. It represents 
a set of two coupled nonlinear kinetic equations for QPs 
and bosons, together with a self-consistency equation for 
the energy gap A. Although, those equations were ac- 
tively studied before 0, H, EM El, H3] , previous studies 
were limited to a simplified version, in which nonequi- 
librium distribution functions of QPs, / = F + Sf, and 
bosons, g — G + Sg, were linearized with respect to equi- 
librium values F and G. As will be shown below, such 
linearization is valid only for extremely small disequilib- 
rium, Sf <C F, Sg <C G, which is hardly satisfied at low 
enough T or high enough energy, where F, G — > 0. 

In this work I present numerical solutions of the full set 
of nonlinear equations, with a focus on strongly nonequi- 
librium phenomena in stacked JJs at high bias. A new 
type of radiative solution is found, resembling lasing in 
semiconducting heterostructures. 

I consider two voltage biased stacked JJs. The role of 
the second JJ is to enhance nonequilibrium QP popula- 
tion by preventing QP leakage from the middle electrode, 
just like in semiconducting double heterostructures [l4[. 
On the other hand, bosons are not confined and can 
freely escape through JJs. The boson escape rate and 
the QP injection rate, relative to the QP relaxation rate 
in electrodes, are described by coefficients "fu and ji, re- 
spectively. All characteristics are shown for the middle 
electrode. Exact formalism, details of the numerical pro- 
cedure and estimation of parameters for Bi-2212 IJJs can 
be found in the Supplementary [13]. 

The QP current in JJs rises stepwise at the sum-gap 
voltage V g — 2A/e. At V > V g it has two sharp maxima 
at E' = and E' = eV — 2A, as shown by the dashed- 
dotted line in Fig. Q]a). Here E' = E — A is counted from 
the edge of the gap. Injected QPs relax via spontaneous 
and stimulated emission of bremsstrahlung bosons and 



FIG. 1: (Color online). QP (a) and boson (b) spectra for 
the ordinary absorptive solution at eV = 3.0Ao. Solid lines 
represent nonequilibrium Sf, Sg, dashed lines-equilibrium F, 
G. The dashed-dotted line in (a) shows the QP injection 
rate. It is seen that there is net accumulation of QPs and 
absorption of bosons at low energies, c) Comparison of QP 
spectra, normalized by the injection rate 7/, for linear (dotted 
lines) and nonlinear solutions for V = 1.75Ao/e < V g and 
V = 2.5A /e > V g and for 71 = 0.1 and 7/ = 0.5 (the red 
dashed-dotted line, for eV = 2.5A only). At eV = 1.75A , 
5/ <JC F and linear and nonlinear solutions coincide. However 
at large bias eV = 2.5Ao, Sf>F and linear solution becomes 
inaccurate. Nonlinearity stimulates QP decay, as seen from 
the strong decrease of 8f(0)/fi with 7/. 



FIG. 2: (Color online), (a) Nonlinear boson emission spectra 
at different V. Bremsstrahlung ^ Q ^ eV — 2A and recom- 
bination £1 ^ 2A bands with sharp maxima at band edges are 
clearly seen for eV = 2.49 and 3.01Ao. (b) Bias dependence 
of Sf(E' — 0) (thin line), A (dashed line, right axis, note 
the reversed scale) and the number of recombination bosons 
5Nbos(Q = 2A) (thick line). Step-like increase of 5Ni, os occurs 
at eV = 4A ~ 3.54Ao upon collision of bremsstrahlung and 
recombination bosonic bands. This also leads to appearance 
of secondary QPs and the secondary bremsstrahlung band 
with < f2 < eV — 4A (downward arrows in a) . 



eventually recombine into Coo p er p airs with emission of 
recombination bosons @, H, EO, [13 ■ Opposite processes 
of boson absorption and pair-breaking slow down the QP 
decay. The rates of all those processes are proportional 
to electron-boson spectral function [13], i.e., QP decay is 
facilitated by electron-boson interaction. 

Fig. Q]represents full nonlinear solutions for a) QP and 
b) boson spectra at V = 3Ao/e, T = 0.5T C and for QP 
injection and boson escape coefficients 7/ = =0.1. In 
Fig. [T] c) linear and nonlinear solutions, normalized by 
7/, arc compared. At low bias V — 1.75A /e < V g the 
QP current is small and 5f <C F. In this case lineariza- 
tion is valid and both solutions coincide. However, at 
V = 2.5Ao/e > V g the QP current increases dramatically 
so that Sf>F and linearization becomes invalid. The 
red dashed-dotted line in c) represents the full solution 
for a larger QP injection rate 7/ = 0.5 at eV = 2.5Ao- 
It is seen that Sf/^/i at the same eV — 2.5A decreases 
substantially with increasing 7/ and disequilibrium, in 
contrast to the linear solution for which Sf/-fj is inde- 
pendent of 77. For 7/ = 0.5, Sf(0)/ji is an order of 
magnitude smaller than for the linear solution, implying 



that nonlinearity tends to stimulate QP decay. 

Fig. [2] a) shows boson emission spectra at differ- 
ent bias. Bremsstrahlung and recombination bands are 
clearly seen at eV = 2.49 and 3.0lAo. Sharp maxima 
at the band edges fl = eV — 2A (upward arrows) and 
il = 2A indicate the most probable QP decay scenario 
Q : first QPs fall to the edge of the gap and then recom- 
bine therefrom into Cooper pairs, emitting = eV — 2A 
and fl — 2 A bosons, respectively. 

Nonequilibrium QPs suppress A from its equilibrium 
value Ao(T). In Fig. [2] a) this is seen as a gradual shift of 
the recombination band edge = 2A with V. Accord- 
ing to the self-consistency equation, the gap is mostly 
affected by QPs at E' = 0. Indeed, from Fig. [2] b) it 
is seen that ^-dependence of A follows <5/(0). Note a 
hysteresis upon sweeping through V g . It is caused by 
nonlinearity of the self-consistency equation. 

Despite substantial quantitative differences, the non- 
linear solution described above is qualitatively similar to 
the linear solution. From comparison of QP injection 
rates Iqp(E) and Sf(E) in Fig. [T]a) it is seen that the 
amount QPs is reduced with respect to Iqp at high en- 
ergies E' ~ eV — 2A as a result of QP relaxation, but is 
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FIG. 3: (Color online). QP (a) and boson (b) spectra for the 
radiative solution. Presentation is similar to Figs. 1 a,b). Re- 
duction of (5/(0) at the edge of the gap in (a) indicates accel- 
erated relaxation due to stimulated boson emission. The cor- 
responding increase of low-frequency boson population 5g(0) 
is indicated in (b). Inset in (b) shows ^-dependence of A: 
gap suppression is very small due to small 8f(0). 



largely enhanced at E 1 ~ 0. Accumulation of QPs at the 
edge of the gap occurs not only because recombination 
is the only sink channel for QPs in the considered case, 
but also because QP relaxation is slowed down by reab- 
sorption of bosons. As a result there is a net absorption 
of bosons at low SI, 5g(0) < 0. Therefore, I will refer to 
this type of solution as the "absorptive" solution. 

With increasing disequilibrium, above some threshold 
QP injection rate, the absorptive solution becomes un- 
stable and switches to a different one, shown in Fig. [3] 
for the same parameters as in Figs. Q]and[U In this case 
there is no significant accumulation of QPs at E' = 0. 
As a consequence, suppression of A is very small, as 
shown in the inset. Instead there is a large net emis- 
sion of low-energy bosons, marked by an upward arrow 
in Fig. [31 b). Apparently, relaxation of QPs is accel- 
erated by stimulated emission of bremsstrahlung bosons 
(note that even for the absorptive solution the nonlinear- 
ity accelerated QP decay, see Fig. Q]c). Significance of 
stimulated emission is obvious from the large total bo- 
son population g(il -> 0) > 1, Therefore, I will refer 
to this new type of solution as the "radiative" solution. 
Once achieved, the radiative solution is sustained even 
at lower bias, typically down to V g , as is the case in Fig. 
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FIG. 4: (Color online). Normalized tunnel conductance vs. 
bias voltage for absorptive (blue) and radiative (red) solu- 
tions and for the same junction parameters. The double-gap 
dip/step occurs in both cases at V = 2V g = 4A/e. How- 
ever, the gap is much stronger suppressed for the absorptive 
solution, despite similar dissipation powers along the I — V 
curves, shown in the inset. Small gap suppression in the ra- 
diative state is due to efficient radiative cooling. 

[3] I want to emphasize that existence of two solutions is 
the consequences of nonlinearity of kinetic equations. 

Upward arrows in Fig. [3] b) and the inset in Fig. 
[3J b) indicate a singularity at the double-gap voltage 
V = 2V g = 4A/e. At this bias bremsstrahlung and 
recombination bands overlap, as seen from spectra at 
eV > 3.54A in Fig. [3 a). This leads to: 

A) Appearance of secondary QPs and a secondary 
bremsstrahlung band: at V > 2V g high energy 
bremsstrahlung bosons become capable of breaking 
Cooper pairs. This leads to appearance of secondary QPs 
with energies ^ E' ^ eV — 4A [9[ . The corresponding 
step-like increase of 5f(0) is clearly seen in Fig. [2] b). 
Relaxation of secondary QPs leads to appearance of sec- 
ondary bremsstrahlung bosons with ^ f2 ^ eV — 4A, 
marked by downward arrows in Fig. [5] a). Recombination 
of secondary QPs leads to enhanced emission of £1 = 2A 
bosons. At eV > 4 A the amount of such bosons increase 
stepwise, as clearly seen in Fig. [5]b). 

B) Double-gap singularity in dl/dV: Step-like increase 
of <5/(0) leads to the corresponding step-like decrease of 
A(V), see Fig. [5]b). This in turn leads to appearance 
of a dip/step in tunneling conductance as demonstrated 
in Fig. |4j Note that the feature is sharper for the ra- 
diative solution due to sharper step in A(V), despite 
smaller absolute gap suppression. Similar features oc- 
cur at all even-gap voltages eV = 2nA, (n = 2,3...) 
due to collision between recombination and secondary 
bremsstrahlung bands 0] . 

Noticeably, I — V curves for both solutions are very 
similar, as shown by the inset in Fig. 2J The large dif- 
ference in A(V) for absorptive and radiative solutions 
for similar dissipative powers P = IV emphasizes in- 
applicability of the concept of thermal conductivity k 
for description of strongly nonequilibrium phenomena in 
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IJJs [2l[. In superconductors n rapidly freezes out at 
T — » 0. Numerical simulations presented here correspond 
to k — because there is no diffusive escape neither for 
QPs nor bosons. Nevertheless, neither QP nor boson sub- 
systems are strongly overheated. This occurs because 
electric power is effectively (with 100% conversion effi- 
ciency) radiated from JJs via ballistic boson emission. 
The phenomenon is well known as radiative cooling in 
light emitting diodes [i~4j . 

I want to emphasize similarities between nonlinear 
nonequilibrium phenomena in stacked JJs and semicon- 
ducting heterostructures [HI], especially Quantum Cas- 
cade Lasers [H, HH- The absorptive and radiative so- 
lutions, obtained here, are similar to light emitting and 
lasing states in injection diodes. Lasing in modern diodes 
occurs above a threshold Jth ^ 10 — 100 A/cm 2 at room 
temperature and Jth decreases exponentially with de- 
creasing T pjj]. For UJSjJjV = V g ) ~ 10 3 -10 4 A/cm 2 at 
much lower T = A.2K [2l|, (23|. According to preliminary 
results (not shown) the threshold for switching from the 
absorptive to the radiative state in JJs also strongly de- 
creases with T. This is caused by the exponential depen- 
dence F (T) which makes it easier to achieve strong non- 
linearity 5f S> F at low T. Furthermore, large Bi-2212 
mesa structures used in emission experiments [l|| HH 
may contain up to 1000 stacked IJJs, which should re- 
sult in significant cascade amplification of boson popula- 
tion (l2l |. Therefore, extreme disequilibrium (population 
inversion) can be achieved in such mesas. 

Obtained results are consistent with several experi- 
mental observations: (i) I have shown that QP relaxation 
becomes nonlinear when 5f > F (see Fig. 1 c). Since 
F freezes out exponentially with decreasing T, QP relax- 
ation inevitably becomes nonlinear at low enough T . In 
this case recombination occurs between two nonequilib- 
rium QPs rather than one nonequilibrium and one equi- 
librium QP, recombination time becomes independent 



of T and dependent of QP injection rate, as observed 
in LTSC JJs at mK temperatures [1, @. (ii) Similar 
dips/steps in dl/dV at even-gap voltages (see Fig. 4) 
were observed for small Bi-2212 IJJs [12j]- Abrupt ap- 
pearance of secondary QPs at those voltages (see Fig. 
2 b) was directly observed in LTSC JJs 0,i,@]. (iii) 
Appearance of radiative cooling upon establishing the 
radiative state in IJJs is consistent with the observed 
small suppression of A and strong reduction of self- 
heating at high bias [2l| . The radiative state can be also 
responsible for recently observed laser-like photon emis- 
sion from large Bi-2212 mesa structures [l5j]. In that case 
the mesa itself acts as a Fabry-Perot resonator, selecting 
cavity (Fiske) modes [tEEl- 

In conclusion, solution of the full nonlinear set of 
kinetic balance equations allowed analysis of arbitrary 
strong nonequilibrium phenomena in Josephson junc- 
tions. This is important for quantitative analysis of su- 
perconducting devices at low T and for understanding 
of nonequilibrium boson emission from IJJs. The discov- 
ered radiative state indicates a possibility of realization of 
a new type of superconducting cascade laser (SCL) [l2l ]. 
Unlike the existing Josephson oscillators which utilize the 
ac- Josephson effect for conversion of electric power into 
radiation [25|, the SCL is based on direct conversion of 
electric power into boson emission via non-equilibrium 
QP relaxation upon sequential tunneling in stacked junc- 
tions. The mechanism is similar to lasing in semicon- 
ducting heterostructures [3, H2] and allows very high 
radiation efficiency. The advantage of employing super- 
conductors instead of semiconductors is in the smaller 
energy gap, which facilitates operation in the important 
THz frequency range. 
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